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SUMMARY 


Integrating matrices allow tha officiant and accurate integration of 
funetiona whose values are given niaaerically on a discrete set of grid 
points. They also form the basis of an efficient numerical procedure for 
solving differential equations associated with the dynasties of rotating 
beams. By expressing the partial differential equations of action in matrix 
notation, utilising the integrating matrix as a spatial operator, and 
applying the boundary conditions, the resulting ordinary differential equa- 
tions can be cast into standard eigenvalue form upon assumption of the usual 
time dependence. As originally developed, the technique has been limited to 
beams having continuous mass and stiffness properties along their lengths. 
This report extends integrating matrix methods to treat the differential 
equations governing the flap, lag, or axial vibrations of rotating beams 
also having concentrated masses. Inclusion of concentrated masses is shown 
to lead to the same kind of standard eigenvalue problem as before, but with 
slightly modified matrices. 


INTRODUCTION 

The equations of motion governing the vibrations and stability of 
rotating beams such as helicopter rotor blades have no closed- form solu- 
tions, and approximate methods of solution such as asymptotic techniques, 
Galerkin's method, or direct numerical integration must be employed. A 
numerical procedure based on the u»e of integrating matrices (refs. 1 and 2) 
has been employed to solve for the vibrations and stability of a wide 
variety of rotating beam configurations with continuous mass distributions 
(refs. 3 and 4). The integrating matrix provides a meats for numerically 
integrating a function that is expressed in terms of the values of the 
function at a set of discrete grid points in the interval ct interest. 

Recent work (ref. 5) has removed the previous restriction thet the grid be 
uniform, and arbitrary increments in the independent variably are now 
permitted, by expressing the equations of motion of the rotating beam in 
matrix notation, utilizing the integrating matrix as an operator, and 
applying the bounds. y conditirns, the resulting ordinary differential 
equations can b<j c'*t into standard eigenvalue form. Solutions can then be 
determined by standard methods. 


As originally formulated (ref. 1) the integrating matrix method 
requires that the mass and stiffness distributions of the beam be at least 
piecewise continuous along the length of the beam. Subsequent applications 
of the method (refs. 2, 3 and 4) were to continuous beams. Hie important 
case in which one or more concentrated masses are located along the beam is 
excluded. Reference 5 allows treatment of discontinuous coefficients in the 
governing equations, but not the case of concentrated masses. As a number * 

of situations of interest involve rotating beams with concentrated masses, a 
need exists to extend the present formulation of the integrating matrix 
technique. 

This report describes generalizations of the integrating matrix method 
required to treat rotating beams with concentrated masses at arbitrary pos i- 
tions along their lengths. Inclusion of concentrated masses is shown to 
lead to a standard eigenvalue problem of the same form as before, but with 
slightly modified matrices. To provide a framework for the generalization 
and to set notation, the following section will briefly review the 
integrating matrix method for vibrations of rotating beams with continuous 
mass and stiffness distributions, 

REVIEW OF THEORY FOR BEAMS WITHOUT 
CONCENTRATED MASSES 

Let ra(x) (0 < x < i ) describe the continuous mass distribution of the 

rotating bean. Consider first tne case of vertical If lap) bending. Let 
w(x,t) be the displacement at position x and time t of the beam normal 
to the plane of rotation. Then, w(x,t) is a solution for 0 < > < i of the 
partial difterential equation (ref. 4): 

[EI(x) w"] M - [T ( x) v' ] f ♦ m(x) v - 0 (1) 

3 a 

where 1 * — , • * — , El is the bending stiffness, ft is the rotation 
3x 3t 

speed, and the tension T is given by: 

* 

i 

T(x) - £l 2 / n m(n) dn (2) 

X 
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Assoc iateu boundary conditions for ths c lamped- free beam srs: 
w(o,t) • w' (o,t) * 0 


(3) 


and 


w"(t ,t) - w" U,t) - 0 (4) 

To oblain the so-called "fundaraental derivative" w", equation (1) may be 
formally integrated twice trom x to t using the boundary conditions in 
equation (4). This gives 

t 

El(x) w"(x,t) ■ - / { m(h ) lw(n,t) - w(x,t)j 

x 

♦ m(n) w(n,t) (n - x)} dn (5) 

Equation (5) can now be converted into an integral equation for the 
fundamental derivative through use ot the relationships 

x 

w'(x,t) * w'(o,t) ♦ / w"(s,t) ds (6) 

o 


and 


x 

w(x,t) “ w(o,t) + / w'(s,t) ds (7) 

o 

As a first step toward expressing equation (5) in matrix eigenvalue 

form, let x , xi , . , . , x^, be N ♦ l grid points along the bean such that 
o N 

U - x < xi <. . . . <. x • i (8) 

o * N 

Spacing between grid points must be non-zero, but need not be uniform. 
Also, let the N + 1-by-N ♦ 1 integrating matrix on this grid be denoted by 
|l] (ret. 5). Then, if { f) is a vector giving the values of a function 
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f(x) at the N + 1 grid points, the action of [L] on {f} can be written in 
the general form 



f (x)dx 


[L] { f } 


(9) 


where i * 0, 1, 2, N. This relation, equations (6) and (7), and the 

boundary conditions of equation (3) now allow w and w' to be re- 
expressed in terms of the fundamental derivative w" through the matrix 
equations 

(w'} * [L] { w"} , 

and 

{w} * [L] 2 (w"} (10) 


where { w} , {w'}, and {w"} are column vectors giving the values of the 
respective functions at the N + 1 grid points x 0 , xj , . .., Xjj. 

To deal with the integral in equation (5), which goes from x to t 
rather than o to x, a slightly modified integrating matrix is required. 
Let [IJ be the N + 1-by-N + 1 identity matrix and let [Bj ] be the N ♦ 1-by- 
N + 1 matrix 


iBi) 



1 


1 


Then, as 


4 l x 

/ f ( x)dx 3 / f(.0dx - / f(x)dx 

x o o 


( 11 ) 


( 12 ) 


the matrix 


[J] » ( t»i I - HD [Ll 


(13) 


4 



ORIGINAL PAGE IS 
OF POOR QUALITY 

is Che required integrating matrix on the grid from x to A, i.e., 


J t(x)dx 


x . 
1 


- U] (f> 


( 14 ) 


where i • 0, . . . , N ♦ 1 . 

Assume the usual tiae dependence 


w(x,t) » w(x) e 


At 


(15) 


o that w' * w' e^ C and w ■ A 2 w e' C . Equation (5) becomes 


- At 


El w" -*■ J m(n) fl 2 n (w(n) - w(x)j dn 


A 2 / m(n ) (x - n ) vtn ) dn 
x 


(16) 


Evaluating equation (lb) at the N ♦ 1 grid points, using the integrating 
matrices lL] and IJ J as operators, and expressing Che resulting equations in 
state variable form leads to the eigenvalue problem 


[G) ) - A i H J {$i ) 


(17) 


where 


{$ 1 ) 


{?'} 


{A w»'}| 


(Id) 


The 2N + 2-by-2N ♦ 2 matrices [ G J and (H] in equation (17) have the 
functional forms 


lG) 


[Gil] 

loj 

[o] 

III 


(19) 
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a^d 


[H] 


Co) 

[Hi2 ) 

II] 

to] 


( 20 ) 


where the N + 1-by-N + 1 matrices IG^J and (H 12 ) are functions of the beam 
properties and the integrating matrix, and [o] is the N 4- 1-by-N ♦ 1 aero 
matrix 

Eigenvalue problems for the in-plane (lag) vibrations and extensional 
(axial) vibrations of a rotating bean are derived in an analogous manner 
(ref, 4). 


INCLUSION OF CONCENTRATED MASSES 


This section describes generalizations of the integrating matrix method 
required to treat locating beams with concentrated masses. The case of 
interior concentrated masses is treated first. The boundary mass situation 
(a tip-mass for the clamped-free beam) is treated by applying a consistent 
limiting procedure to the interior mass case. Flap, lag, and axial vibra- 
tions are treated separately for convenience of discussion. 

Interior Masses. 

Consider the effect of a single concentrated mass ot magnitude M 
located at position o\5\£on the flap, lag, and axial vibrations 

ot a rotating beam. This entails no loss ot generality since the eftects of 
multiple concentrated masses may be commuted separa' 'Jy and then sunraed to 
obtain an aggregate ettect. In what follows, the total mass distribution of 
the beam (distributed plus concentrated) will be denoted by U(x). Thus, if 
6(x) is the usual Dirac delta function and m(x) again denotes the dis- 
tributed mass 

u(x) ■ ra( x) ♦ M <5(x - £ ) . (21) 
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Vertical (Flap) Vibrations . The motion of the rotating beam with con- 
centrated inboard mats is still governed by equation (1) if m is replaced 
by the total mass distribution u . Further, as o < £ < t, the boundary 
conditions of equations (3) and (4) still hold, thus, the derivation of 
equation (5) remains valid if m is replaced by y. The Dirac c*elta 
function has the property that 

1 f(£) if x < t < i 

/ 6(n - O f(n) dn » (22) 

x o if o < 5 < x 


Equation (5) now becomes 


£ 

El w ,, (x l t) * -/ { m(n ) fl 2 n lv(n,t) - w(x,t)] 
x 

♦ m(n) w(n,t) (n - x)} dn 

- D^(x) {m lw(C »t) - w( x , c ) | r M vU ,c)(C - x)} 


(2 3) 


where 


l\(x) 


U if o < ( < x 
1 if x < ( < t 


(24) 


Assuming a solution of the form v(x,t) * w( x) e gives as the analogue of 
equation (lo) the equation 

£ _ 

El w 1 ' ♦ / m(^ ) ft 2 n Iw'n ) - w(x)J dn + Mft 2 t D (x) (w(£) * w( x ) ) 
x ** 

£ 

* A 2 {/ m(n)(x - n) w(n ) dn ♦ MD (x)(x - t) w(£)} (25) 

x ^ 


The effect of the concentrated mass M at x • £ enters equation (25) only 
through the additive terms 


(x U (x) («(*> - w( x ) j 


(26) 
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X 2 H^x.c) - x 2 MD (x)(x • e) v XO 


(2 7) 


with the remainder of equation (25) being exactly the same as in the 
previous case of distributed mass only. 

To express equation (25) in matrix form, again consider the discrete 
grid of N + 1 points 

0 ■ x < xi < . . . < x k * l 
o 1 N 


where £ need not coincide with a grid point, and let [L] be the 
integrating matrix on this grid. Consequently, when equation (25) is 
evaluated at the grid points and expressed in matrix form, the resulting 
eigenvalue problem in state variables will have the form of equation (17) 
with modified matrices 


[G] 


tG U ] + [Gill 
[ol 


[0] 

[1] 


and 



[Hi 2] + (Hi 2) 

to] 


(28) 


(29) 


where [Gi i ] and [Hi 2 1 are the same matrices as in relations (19) and (20), 
and [o^] ard [Hjji are N + 1-by-h + 1 matrices which gi /e the values of 
Gj (x,£ ) and Hj (x,£ ) at the N + 1 grid points. 

The principal difficulty in deriving 1 J and Ud l2 J is that both 
Gj (x,£ ) and Hj (x,£ ) involve w(£) but 5 need not be a grid point. 

However, this problem may be overcome by using an interpolation formula 
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(such aa Lagrange interpolation) on the grid with coefficient a that depend 
only on £ and the grid pointa, x 0 , . .., x(j. In particular, the 
coefficienta ahould not depend on valuea of the function to be interpolated 
at the grid pointa. The interpolating expreaaion will thua be of the form 


wU) ■ * 0 (5) w(x q ) ♦ wixj) + ... ♦ wix^) (30) 


For example, if an Mth order Lagrange interpolating polynomial on the aubaet 

of grid points x , M i* chosen, where x < x < 5 

Y Y + 1 Y + M o Y 

< x^ u < x , then in equation (30) 

T + M N 

a, (5 ) = 0 for k < Y or k > Y + M 
k 


while for Y * k < Y+M 

s (t , . ••• «- Vl ) «-Vi> ... 

k <vV ••• 'VVi 1 'VVi 1 'vs * 1 


( 31 ) 


By choosing M sufficiently large, e.g. , M ■ 7, the resulting interpolating 
polynomial gives a high degree of accuracy without the need to cluster grid 
points around £. To use equation (30) in the present context, let IA^ J 
be the N + l-by~N + l matrix with ik-th element 

(A_K - a ,(S) k - 1, ..., N + 1 (32) 

5 lk k-1 


i.e., [A^ ] has identical rows and each row has as elements the successive 
coefficients in the interpolating formula (30). Pre-mul tiplying (w) by 
[A^J, equation (30) now immediately shows 

[A ] - { w(5 )) (33) 


Some care must be taken in defining a matrix that gives the effect at 
the grid points of multiplying by the discontinuous function Dg(x) which 
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change* fro* unity to aero aero** * ■ I. Suppose that C lie* in the 
subinterval 


x j - x j • L, ...» N 


( 34 ) 


Then, if 
entries 


i* the N + 1-by-fl + 1 diagonal matrix with diagonal 



(35) 


(3b) 


then 


[V] tf) *t 


gives the vector of values ,of the function D$(x) f(x) 


H 


{ f} does not 


at the grid points. However, it should be noted that 

explicitly reflect the non-aero values of D^(x) f(x) on” the interval x. , 

£ J - 1 

< x < Consequently, in specific applications, the grid should be chosen 
if possible so that £ - xj « 1 is small relative to xj - £ , 


Let the N ♦ 1-by-N ♦ 1 diagonal matrices 


by 


and N be defined 


N- 


II] 


(37) 


and 




(i - 1 N ♦ 1) 


(3d) 
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Equation* (26) and (27) now (mediately show that the matrices (Cul 'nd 
[Hj 2 1 have the forms 


(Cll] 


[\] |\] t'V - ll|] < L|S 


(39) 


and 

I?i2l • [\J [\J IV IL1 2 


(4U) 


In-plane (lag) vibrations . Lee v(x,t) denote the deflection of the 
beam in the plane of rotation, where the total mass distribution is again 
given by y(x) in equation (21). Seeking solutions of the form 

v(x,t) • v(x) e^ C now gives the governing differential equation 


El v"(x) ♦ / m(n) n 2 {nlv(n) - v(x)J - 7(n) In - x] > dn 


+ d^(x) m 2 {£[7(0 - v(x)j - 7(£) (£ - >ij > 


- X 2 


It 

/ m(n ) lx - n ] v(n ) dn 

X 


( 41 ) 


+ (x) M(x - £ ) v(£ ) 


With the exception of the additive contributions 


1*2 (x,£ ) ■ *•»! (x,£ ) - L^(x) tt} 2 v(£) l£ - xj 


(42) 


and X 2 Hj(x,£), where Gi(x,£) and Hi(x,$) are defined by equations (2b) 
and (2 7) with w replaced by v, equation (41) is identical to the 
corresponding equation without concentrated masses. An arguaent similar to 
the flap case now gives that the eigenvalue problem associated with (41) 


for X and 


(v") 

{Xv") 


is of the same form as (17), (28), and (29), with 


(Gj j ] and (H 12 ) as in the lag problem without concentrated mass, (H|2 ) again 
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given by ( 40 ), and [G^] given by Che expression in ( 39 ) plus the additional 
tern 


* M M *v il|2 


( 43 ) 


Extensional (Axial) Vibration* . Let u(x,t) be the axial extension of the 
beam and seek solutions of the form u(x,t) * ut( v ' e^ C . Then, the governing 
differential equation for the rotating beam with concentrated interior mass 
is 


l 

AE u ’ - J m(n)G 2 «»(n) dn “ D. (x) rtl Z u(0 

X * 


• -A 2 


l 

I 

x 


m(n) u(n ) dn ♦ d^(x)m u(C) 


I 

f 


(44) 


This equation differs from the case without concentrated mass only through 
the additive terms 


G3 (x, , ) - -D^ (x) Kl 2 u(£) 

and 


(45) 


A 2 H 3 (x,C> - -X 2 (x) M u(r) 


(4b) 


Consequently, the corresponding eigenvalue problem for 


X and ($3) 


Cu’} 

X{u*} 


will be of the same form as (17), (2D), and (29) where [Gnl and [H11) are 
as in the axial problem without concentrated mass, and, as (u) * [L] {u 1 }. 


[Hi 2 ) - -M 



[L], 


(47) 


and 


[Gi 2 1 ■ 0 2 [Hi 2 I • 


(48) 
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Boundary Masses 

If Che concentrated mass M is placed directly at the free end of a 
rotating cantilever beam, i.e. 5 “ t, one of the free-end boundary 
conditions is replaced by th~ equation of motion of the concentrated mass. 
For example, in the flap case the boundary conditions at x ■ l become 

w” - 0 (49a) 


and 


(El w”)' - M(n 2 * w' + w) - 0 V49b) 

An appropriate limiting procedure may now be used to show that equation (23) 
for the interior case yields the conditions of equations (49) as % tends 
to A . In this connection, let Lim denote the limit as x -► A with £ > x. 

This limiting procedure thus automatically ensures that £ ♦ A as x + A . 
Further, as £ > x, (x) = 1. 

Applying the limiting procedure to equation (23) immediately gives 
Lira (El w 11 ( x , t ) ) “ El w M (A,t) * 0 

consistent with (49a). Differentiating (23) with respect to x with £ > x 
gives 

/ , * 

( El w"(x, t))' * - / { m n w'(x f t) - m w(n , t)} dn 

x 

+ M w r ( x , t ) ♦ w(C,t)} (50) 

Applying Lim to (50) now gives (49b). This shows that equation (23) remains 
valid as £ ♦ A . This, in turn, implies that the integrating matrix formula- 
tion for flap vihrations with interior concentrated mass remains valid in 
the case of a boundary mass. Similar arguments, which will be omitted, 
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involving equation (41), (44), and the limiting procedure shot# that inboard 
concentrated mass formulations for the lag and axial vibrations also hold in 
the boundary mass case. 

VERIFICATION OF THE FORMULATION WITH 
CONCENTRATED MASSES 

To verify the present formulation of the integrating matrix method 
including the effects of concentrated masses, a nunber of beam problems were 
analysed. These included the flap, lag, and axial vibrations of beams with 
a single inboard concentrated mass, a tip mass only, and up to ten concen- 
trated masses pi seed along the length of the beam. Both rotating and non- 
rotating situations were considered. Results for these test cases were 
compared to known exact solutions and solutions obtained using a finite 
element program. In all cases, results obtained from the integrating matrix 
formulation were found to be in excellent agreement with these other 
results. For the flap and lag vibrations, test cases also included beams 
with ten concentrated masses of appropriate size and no distributed masses, 
i.e., m(x) i 0. Excellent agreement with other results was also obtained 
for these case9. It is worth noting that the integrating matrix formulation 
is not appropriate tor the axial vibration of a beam with only concentrated 
masses if u 1 is taken as the dependent variable as is done here. Indeed, 
if m(x) « 0, equation (44) has the functional form AE TP - constant, and the 
eigenvalue character of the problem l$ lost. 
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